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Abstract 

We present a study of the B spectrum performed in the framework of Heavy Quark Effective 
Theory expanded to next-to-leading order in 1/mb and non-perturbative in the strong 
coupling. Our analyses have been performed on Af = 2 lattice gauge field ensembles 
corresponding to three different lattice spacings and a wide range of pion masses. We 
obtain the Bg-meson mass and hyperfine splittings of the B- and Bg-mesons that are in 
good agreement with the experimental values and examine the mass difference — m-Q 
as a further cross-check of our previous estimate of the b-quark mass. We also report on 
the mass splitting between the first excited state and the ground state in the B and Bg 
systems. 


Keywords: Lattice QCD, Heavy Quark Effective Theory, Bottom mesons, Spectroscopy 
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1 Introduction 


Lattice studies of the hadron spectrum provide valuable indications on the reliability of ef¬ 
fective held theories. A phenomenologically relevant example is given by effective held the¬ 
ories for heavy quarks, such as Heavy Quark Effective Theory (HQET) or Non-Relativistic 
QCD (NRQCD). In this respect, hyperhne splittings are of great theoretical interest, since 
they vanish in the static (1/mh — >• oo) limit, and thus serve to probe subleading terms in 
the heavy quark expansion. 

Regarding specihcally the hyperhne splitting of the B spectrum, it is well known 
that the masses of the pseudoscalar Bg and vector B* differ as a result of spin effects [1]. 
This mass difference is produced in HQET by a 0(l/mh) chromomagnetic term in the 
Lagrangian that breaks the heavy quark spin symmetry of the static theory [2,3]. In the 
context of lattice NRQCD it was found that in order to obtain a determination of the 
hyperhne splittings of bottomonium that is consistent with experiments, the perturbative 
improvement of the NRQCD action to at least one-loop order is required [4-6] , for values of 
the lattice spacing of about 0.1 fm. The T spectrum, on the other hand, is well reproduced 
within NRQCD, and has been used in the literature as a quantity to set the lattice spacing 
[7,8]. Hyperhne splittings have also been studied within the HISQ approach to heavy 
quarks [9]. 

In contrast to other studies, our results are based on non-perturbatively renormalized 
HQET including the next-to-leading order terms in the inverse heavy quark mass expan¬ 
sion. When formulated with a lattice regulator, the resulting lattice held theory can be 
non-perturbatively matched to QCD in the continuum limit. This provides a rigorous and 
systematically improvable approach to the study of the B-meson system. 

Although permil precision seems difficult to reach within this approach, as it would re¬ 
quire including 1/m^ corrections, the control over systematics such as the non-perturbative 
subtraction of power-divergences and the related existence of the continuum limit, makes 
the approach very appealing and conceptually sound and in this sense, and in our view, 
more precise than other methods. In addition, as evident in [10] and emphasized in [11], 
for quantities more “complicated” than decay constants, such as form-factors, there is great 
need for results from different approaches, with different levels of control on the various 
systematics. Applications of the non-perturbative HQET formalism to form-factors are 
well on their way [12-14]. 

Remaining within spectrum-related observables, another quantity of interest is the 
SU(3) isospin breaking difference msg — tu-q. Most of the dependence on the heavy quark 
mass cancels in the difference and the statistical correlations between the strange and 
light-quark measurements should lead to a much improved precision in the determination 
of the difference compared to what would be possible for the individual masses. In [15], the 
mass of the B-meson has been used to determine the b-quark mass. Since the mass of the 
Bg-meson could have equally well been used, it is instructive to study how the strange-light 
mass difference propagates to the b-quark mass measurement. 

Relatively little is known experimentally on the radial excitations of B-mesons, usually 
denoted by B'. CDF has claimed the observation of a resonant state B(5970) that may be 
interpreted as a B' state [16]. Lattice predictions for the radial excitation energies of the B 
system have been made in the framework of NRQCD [17,18] and HQET [19-27], but, with 
the exception of the most recent HQET studies, control of the continuum extrapolation has 
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/3 = 

5.5 

/? = 

5.3 

/3 = 

5.2 


HYPl 

HYP2 

HYPl 

HYP2 

HYPl 

HYP2 

®^bare 

0.969(12) 

1.000(12) 

1.317(15) 

1.350(15) 

1.520(17) 

1.553(17) 

^^bare 

0.594(14) 

0.606(14) 

0.993(17) 

1.014(17) 

1.214(18) 

1.239(18) 


0.520(11) 

0.525(10) 

0.415(8) 

0.418(8) 

0.377(7) 

0.380(7) 

i^spin/^ 

0.949(37) 

1.090(42) 

0.73(28) 

0.883(33) 

0.655(24) 

0.812(30) 


Table 1: Values of HQET parameters at the physical point uj(z = Zh)- As determined in [15], we 
have used = 13.24 to interpolate the parameters of HQET at static order, and Zb = 13.25 for 
the parameters of HQET expanded to NLO. The bare coupling go is given by /3 = B/gg. 

been limited, if possible at all. An additional difficulty that always needs to be addressed is 
in disentangling single particle excitations from multi-particle states. We use the notation 
B'" (or with g = u/d or s) for the excited states that we observe. They might be 
identified with radial excitations B', but without a dedicated study including also multi¬ 
hadron operators, we are unable to conclusively determine the nature of these excited 
states. 

In this paper, we report on our estimate of the hyperfine splitting mp* — tub in the B 
and Bg systems, and on the mass differences niB^ — ut-b, w-bj ~ and thb^ — tub from 
Nf = 2 lattice simulations. 


2 Theoretical setup 

2.1 HQET on the lattice 

The HQET action at 0(l/mh) reads 

'S'hQET — n ^ (^) ^kin^kinix') ^.aspinOspin(^)} j (^T) 

= 'iphix) Do iph{x), (2.2) 

C’kin(a^) = V’h(a:) V’h(a;), (2.3) 

C’spin(a;) = 'iphix) o- ■ B 'i/jh{x) , (2.4) 


where the subscript h denotes a heavy (static) quark field satisfying = V’hj and the 

operators are normalized such that the classical (tree level) values of the coefficients are 
<^kin = Wspin = l/(2mh). The energy levels computed in this theory are relative to a bare 
mass mbarej which at tree level is simply m-h, but has to absorb a power-divergent shift 
at the quantum level, implying that it will take a different value depending on whether 
the 0(l/mh) terms are included or not. The renormalizability of the static theory can be 
preserved at NLO by considering a strict expansion of correlation functions in 1/mh, 

{O) — (O)stat T Wbin O ^ (O 0kin(^))stat T ^spin ^ ^ ^ ^ (O Ospin(a;))stat ) (2-5) 

where the suffix “stat” denotes an expectation value measured in the static theory. 

The values of the parameters and {mbare; ^^kim <^spin} required to match to QCD 

have been determined in [15,28] using the Schrodinger Functional scheme. In order to keep 
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this paper self-contained, we summarize in Table 1 their values for the three lattice spacings 
o(/3) and two static discretizations (HYPl, HYP2 [29]) that we use. The parameters are 
given at = TiM^ where Mb is the Renormalization Group Invariant (RGI) b-quark 
mass and Li is the linear extent of the lattice used for the matching. 


2.2 The variational method 

For spectroscopic applications, having good control over contributions from excited states 
is crucial to reduce the systematic error on the determination of the ground state. For this 
purpose it is necessary to use a variational method [30-34], which starts from matrices of 
correlation functions, 

cr(i) = E + i,y)OK^))stat> 

(2.6) 

^kin/spm^^^ = E (^*(^0 + t, y) O* (x) C’kin/spin(^)),t,,t , 

containing all pairwise correlators of some set of interpolating fields Oj, i = 1,..., N, and 
proceeds to solve the generalized eigenvalue problem (GEVP) 

^stat(^) ^stat(^^ ^ to), n = 1, . . . , iV , t > to . (2.7) 


From the generalized eigenvalues to) and the corresponding eigenvectors v^^^{t, to), 

the effective energy levels can be computed as [34,35] 


Ef’^^-\t,to) 

to) 

where x G {kin, spin} and 


-^[logA^*'^*(t + a,to)-logA^*'^*(t,to)], 

Kit, to) _ X^{t + a,to) 

X^^-^it,to) X^^-^{t + a,to)' 


[|Ar(«, *o)i-‘c“(t) -c‘((„)K“(i.(„)). 

The corresponding asymptotic behaviour is known to be 

Ef’^^^\t, to) = E^^^^ + ' + ..., 

EfKt,to) = E^ + [K- /3r + , 


( 2 . 8 ) 


(2.9) 


( 2 . 10 ) 

( 2 . 11 ) 


where the energy gap is defined as AEm,n = E^ — En¬ 
in our study, the operator basis used is given by heavy-light bilinears with varying 
levels of Gaussian smearing [36] applied to the light quark field, 

Ofc(x) = V'h(a^)7075V’f^(a;), (x) = {l + kg K 'iljq{x), (2.12) 

where g = u/d or s. The covariant Laplacian A is built from gauge links that have been 
triply APE smeared [37,38] in the spatial directions, and the smearing parameters kq and 
Rk are chosen so as to approximately cover the same sequence of physical radii at each 
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Figure 1: Illustration of typical plateaux for the 0(l/77ib) chromomagnetic energy (top) 

and the first excited state static energy (bottom), computed on the CLS ensemble N6 

(a = 0.048 fm, = 340 MeV). 


value of the lattice spacing, as discussed in [15]. We solve the GEVP for the matrix of 
correlators in the static limit for N = 3. 

An illustration of two typical plateaux of the energies and is shown in 

Fig. 1. The plateau regions used for weighted averaging have been chosen by applying the 
procedure already discussed in [39,40] in order to ensure that the systematic errors due to 
excited-state contributions are less than a given fraction (typically 1/3) of the statistical 
errors on the GEVP results. As a consistency check, we have also employed a global fit of 
the form of eqs. (2.10) and (2.11) to our data. The values of En obtained from the fit are 
compatible with the plateau values, and generally exhibit smaller statistical errors. We 
therefore consider our errors to be estimated conservatively. 

2.3 Computing masses 

We start by recalling the definition of the B,j-meson mass as it has been used to non- 
perturbatively fix the HQET parameters for the static quark discretizations HYPl and 
HYP2 in [15]. On the lattice one combines the HQET parameters (see Table 1, or [41] for 
a more detailed discussion) with the large volume computations of ground state energies, 

mgf(m.,a) =mtfe + ^rL,, (2.13a) 

me, (rrin, a) = mbare + Elf'"* | + WkinE^i “ | + uJspinEl^"^ I . (2.13b) 
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We use the subscript mq to indicate that the corresponding energy results from the GEVP 
analysis of correlation functions in the heavy-light {q = d) or heavy-strange {q = s) sector 
along the lines presented in Section 2.2. The PDG value me = 5.2795 GeV in [1] has already 
been used as input to determine mb [15] and thus the HQET parameters Wj. However, re¬ 
computing ms at fixed values of the HQET parameters serves as a non-trivial cross-check 
of our calculation. Additionally, we will compute hr the very same way. All remaining 
observables presented in the present paper are mass splittings which can be computed with 
the data resulting from our GEVP analysis. We now give their explicit dehnition before 
discussing any combined chiral and continuum extrapolation in Section 2.4. 

The hyperhne splitting of the Bg-meson system, me* — w-b^, is given by 

A™m(m^,a) =, 9 = s,d, (2.14) 

and depends on the HQET parameter Wspin and the large-volume measurement of 
The strange-light mass difference — niB^, at the static and 0(l/mb) level, reads, 

(m^, o) = As-dEf ^^, (2.15a) 

As_dm (m^, o) = As-dEf^^ -k WkinAg-dE^i” + WspinAs-dE^i*’” , (2.15b) 

respectively, where we have dehned the shorthand Ag-dQ = Q\mq=mB ~ Q\mq=md ' Finally, 
the mass gaps for excited states, mB^ — mB,, can be computed from AEm,n = E'm ~ En 
via 


A®“m(m^,a) 

Ai.,qm(m^,a) 


AE;|*f I 
AElff I 


+ <^kin‘i^-E'2,T I m ‘^spinAPl^p^ |^ 


(2.16a) 

(2.16b) 


to static and 0(l/mb) accuracy, respectively. As any observable computed on the lattice, 
these quantities depend on the input parameters used in the simulation. The strange 
and bottom quark mass have been hxed to their physical value along the lines reported 
in [15,41,42], such that only the dependences on the light quark mass—here parameterized 
by —and lattice spacing a remain. 

We have extracted the energies on a subset of gauge ensembles from the Goordinated 
Lattice Simulations (GLS) effort for Nf = 2. The parameters of the simulations and 
the statistics entering our analysis are collected in Table 1 of [15]. Remaining details 
of the measurements have been summarized in Table 2} The light quark is treated in 
a unitary setup with 171^ in the range [190 MeV, 440 MeV], and the bare strange quark 
mass has been tuned on each GLS ensemble to its physical value by using the kaon decay 
constant /k to set the scale [42] and tuk = 494.2 MeV. The lattice spacings are o/fm G 
{0.048,0.065,0.075} for /3 G {5.5, 5.3, 5.2}, corresponding to the GLS ensemble ids N-0, 
E-G and A-B, respectively. All lattices have mT^L > 4, such that hnite-volume effects are 
sufficiently exponentially suppressed at the level of accuracy we are working at. 

Our data analysis takes into account correlations between different observables as well 
as intrinsic autocorrelations of the Hybrid Monte Garlo (HMG) algorithm resulting from 
slow modes in the simulation. As these contributions rapidly grow towards the continuum 

^ For the B6-ensemble a preliminary value of Kg has been used, which produces a bare subtracted quark 
mass which differs by about 1 MeV from the Hnal one. 
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e-id 


Kd — /^sea 

Ks 

T 

Tcfg 

^exp 

heavy-light 

HYPl HYP2 

heavy-strange 

HYPl HYP2 

A4 

5.2 

0.135900 

0.1352850 

2 

8 

106 

1012 [1,1] 

1000 [1,1] 

A5c 


0.135940 

0.1352570 

2 

4 

39 

501 [1,1] 

500 [1,1] 

A5d 

B6 


0.135970 

0.1352390 

2 

2 

39 

500 [1,1] 

636 [1,1] 

636 [1,1] 

E5 

5.3 

0.136250 

0.1357770 

4 

16 

151 

1000 [1,1] 

1000 [1,1] 

F6 


0.136350 

0.1357410 

2 

8 

151 

500 [1,1] 

600 [1,1] 

F7 


0.136380 

0.1357300 

2 

8 

151 

602 [100,1] 

200 [1,2] 

G8 


0.136417 

0.1357050 

2 

2 

56 

410 [1,1] 

- 

N5 

5.5 

0.136600 

0.1362620 

0.5 

8 

488 

477 [300,1] 

- 

N6 


0.136670 

0.1362500 

2 

4 

108 

950 [2,2] 

707 [100,2] 

07 


0.136710 

0.1362430 

2 

4 

108 

980 [1,1] 

490 [1,2] 


Table 2: Measurement details for observables in the light- and strange-quark sector. For each 
ensemble we list the light- and strange-quark hopping parameter, the trajectory length r, the 
distance of saved configurations Tcfg and the exponential autocorrelation time Texp in MDU. For 
each individual static quark discretization we then specify N^s i.e., the total number of 

measurements starting from the to**' saved configuration using every d**'. Note that for A5 we 
have two independent Monte-Carlo chains. 


limit, it is mandatory to estimate and include a-priori unknown long-tail contributions 
from the autocorrelation function of each individual observable. Further details of our 
analysis method can be found in Appendix A. 


2.4 Extrapolation to the physical point 


To extrapolate our data to the continuum limit and to the physical point, we take ex¬ 
pressions from heavy meson chiral perturbation theory (HMyPT) if available, and a linear 
ansatz otherwise. 

In the chiral regime, the mass of the B-meson can be extrapolated to the physical 
point using a functional form motivated by HMyPT [43]. Defining a subtracted mass by 
removing the leading non-analytic (in the quark mass) term, viz. 


«) = a) -k c, 


16vr V /, 


m:; 


(/rp)2 ) 


(2.17) 


with Cq = 1 in HMyPT at NLO for (7 = d but zero otherwise, the parameterization reads 

a) = Bq + Cq {y - y®^P) -k . (2.18) 

In (2.17) the B^Byr-coupling 'g = 0.492(29) is taken from [44] and y = m^/( 87 r^/^) (with 
the convention = 130.4 MeV and = 134.98 MeV). We use the same set of 

measurements for /j^ and m.,^ on each CLS ensemble as reported in foregoing analyses [15, 
41]. We add the subscript 5 to distinguish the two available static discretizations which 
are combined in the extrapolation to obtain the parameter Bq = tob, at the physical point 
(y, a) = (y^'^PjO). For the Bq-meson mass an extrapolation quadratically in the lattice 
spacing is justified as we have full 0(a) improvement at work in (2.13) at the static order 
and 0(a) terms are then suppressed by a factor 1/mb once NLO terms in HQET are taken 
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into account. For the case of me that has been explicitly checked in [15] and it is therefore 
conceivably valid also here for mBg. 

In Fig. 2 we show the continuum and chiral extrapolations of and m-Q^ — m-Q at 
next-to-leading order in l/rub. In this and the following Figures, hlled symbols and dashed 
curves represent our HYPl data set, while open symbols and dash-dotted curves represent 
our HYP2 data set. For both, equal colours and symbols refer to the same lattice spacing 
as indicated by equal values of the bare gauge coupling, given by /3 = G/f/g in the legend, 
c.f. Tables 1 and 2. The solid black line is the continuum limit given by the part of the ht 
function which is independent of the lattice spacing a. Accordingly, the difference between 
the solid and non-solid lines represents the cutoff effect as htted through a given ansatz. 

For the hyperhne splitting, we can perform the chiral and continuum extrapolation 
using 

A^Jm {y, a) = [mB* - rriB,] [l - _ ^exp j^^exp) ^c,iy- y-P)] 

+ Dq S a Dq S 0 ? , (2.19) 

with the continuum part coinciding, in the case (? = d, with the expression derived in [45]. 
Hence, we can probe two ansatze for the chiral extrapolation by setting cj = 0,1. Since 
in principle the 0(a) improvement of the hyperhne splitting is not implemented, we have 
included linear cutoff effects, but also study the scaling behaviour with an O(a^) ansatz by 
setting either Dg s or 1)^ ,5 to zero in the equation above. In general, our data is not sensitive 
enough to clearly separate 0(a) scaling from O(a^) such that individual hts lead to similar 
results. As additional safety measure we account for a systematic error by increasing the 
uncertainty of the favoured 0(a) extrapolation in order to cover the mean value obtained 
from the corresponding O(a^) extrapolation. As will be seen in the following, this only 
occurs in the case of ms* — rriBs ■ 

Since, to our knowledge, no systematic HM^PT formulae exist in the literature for 
the mass splittings tub^ — we again employ a simple linear ansatz in m^: 

(y, a) = [mn^ - + Cf'^*(y - y^P) + E^fa^ , (2.20) 

Ar,gm (y, a) = [m-B^ - tubJ + C'(y - y®^P) -F Eg ^a + ^g,sa^ ■ (2.21) 

For the mass splitting A^^gm at next-to-leading order in HQET, one can in fact dehne 

(?/>«) = (y, a) - (y, a) (2.22) 

and take the continuum limit of the leading (static) and next-to-leading order piece (1/m) 
separately. While for the former eq. (2.20) can be employed, the 1/mb-correction term has 
to be extrapolated using ansatz (2.21) with E”^ = 0. Both pieces can again be combined 
after the physical point is reached, leading to an improved signal that is dominated by 
the static extrapolation. No matter what way is pursued, all extrapolations lead to results 
that are well compatible within errors. 

We perform the full analysis using two static discretizations (HYP1,HYP2) while only 
one may already suffice. The reason is that the universality of the continuum limit, a 
property that can be rigorously guaranteed for renormalizable theories, implies that both 
static discretizations have to lead to equal results in the continuum. At hnite lattice 
spacing the two discretizations have inherently different lattice artifacts, best seen among 
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Figure 2: Chiral and continuum extrapolation of mB^ (left) and — tob (right) according 
to (2.18). The open triangle represents the corresponding result for extrapolating the static order 
data. Here and in the following Figures, filled symbols and dashed curves represent our HYPl data 
set, while open symbols and dash-dotted curves represent our HYP2 data set. The solid black line 
is the continuum limit for the given fit ansatz. 

our observables in Fig. 3. To this end one can either take the two independently in order 
to show that universality holds or—as we usually do—combine both in one fit in order to 
stabilise the x^-mhrimization. It furthermore serves as a non-trivial check that the whole 
procedure has been applied correctly and consistently. 

3 Results 

Using the HQET parameters of Table 1 at Zb = 13.25, we obtain for the B-meson mass 
rriB = 5.285(62) GeV by employing eqs. (2.17) and (2.18) with Cg = 1 as extrapolation 
ansatz. Note, however, that the experimental value mp = 5.2795 GeV has been used as 
input in [15] to fix the b-quark mass. Therefore mp is not a prediction of the theory in our 
setup, and the number above should be regarded as a consistency check of the approach. 

3.1 Ground states 

As first quantities beyond tub, we compute the Bg-meson mass and the mass difference 
rriBg — m-B ■ Their continuum and chiral extrapolations at next-to-leading order in 1/mb are 
shown in Fig. 2 and compared to the extrapolated value of the static data at the physical 
point. The raw data for the mass difference mB^ — ut,b and the hyperfine splittings is 
collected in Tables 3 and 4. 

From the extrapolation of the form of eqs. (2.17) and (2.18) with Cq = 0, we obtain 
the result: 

mB, = 5383(63) MeV , mB, - mgf = 58(12) MeV , (3.1) 

where the error includes statistic and systematic uncertainties (scale setting, HQET param¬ 
eters), combined in quadrature as explained in Appendix A. Although our result for rriBg 
comes with a much larger error compared to the PDG value, mBO = 5366.77 ± 0.24 MeV, 
the difference between the mean values is only one fourth of our error. 
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Figure 3: Chiral and continuum extrapolation of the hyperfine spin splitting of the B-meson (left) 
and of the Bg-meson (right). For tob* — ruB we show the extrapolation from ansatz (2.19) with 
Cd = 1 and leading a-effects, while for mB* — mB^ we have to set Cg = 0. In both panels the 
open triangle at the physical point reflects the corresponding fit result using an a^-scaling ansatz. 
Details are given in Table 3. 




Ag dm [MeV] 


[MeV] 

e-id 

y 

HYPl 

HYP2 

HYPl 

HYP2 

A4 

0.0771(14) 

55.2(2.6) 

56.0(2.3) 

51.9(2.1) 

53.1(1.8) 

A5 

0.0624(13) 

68.8(1.8) 

66.0(1.6) 

62.6(1.8) 

63.4(1.7) 

B6 

0.0484(9) 

71.6(4.9) 

70.2(3.5) 

66.5(3.9) 

67.3(3.1) 

E5 

0.0926(15) 

52.8(2.7) 

51.9(1.9) 

48.8(2.0) 

49.2(1.7) 

F6 

0.0562(9) 

70.9(5.7) 

70.3(4.0) 

66.5(5.2) 

64.7(3.9) 

F7 

0.0449(8) 

69.7(3.2) 

74.0(2.5) 

67.1(3.0) 

69.4(2.4) 

N6 

0.0662(10) 

64.0(6.7) 

68.0(3.7) 

62.8(4.8) 

63.3(3.2) 

07 

0.0447(7) 

- 

72.2(6.9) 

- 

70.2(5.3) 

LO-a2 

j/®’^P,a = 0 

91.2(5.8) 

87.2(5.2) 

NLO-a^ 

= 0 

88.9(5.7) 

85.0(5.1) 


Table 3: Mass splitting between Bg- and Bj-meson at static and next-to-leading order HQET. 


In the combination mBg — m-B, the error is reduced due to correlations among the 
heavy-light and heavy-strange measurements, although the latter have been performed on 
a subset of the available ensembles only, as it can be inferred from Table 2. 

Our results for the Bg-B mass splitting at 0(l/m) and in the static approximation are 


mB, -mB = 88.9(5.7)(2.3)x MeV , (3.2a) 

[niB, - mB]"*'"* = 85.0(5.1)(2.2)^ MeV , (3.2b) 

both in good agreement with the PDG value of [1] niB, — tub = 87.35(23) MeV. Here, 
the quoted mean and statistical error results from an HM^PT extrapolation ansatz that 
is obtained by appropriately combining eq. (2.18) for the light and strange quark sector. 
As systematic error due to the chiral extrapolation ansatz we quote its difference to the 
standard m^-extrapolation. 
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[MeV] 

Af^m [MeV] 

e-id 

y 

HYPl 

HYP2 

HYPl 

HYP2 

A4 

0.0771(14) 

52.3(2.1) 

62.0(2.4) 

54.2(2.2) 

63.1(2.5) 

A5 

0.0624(13) 

51.1(2.0) 

60.5(2.4) 

52.9(2.1) 

62.9(2.5) 

B6 

0.0484(9) 

50.1(2.2) 

58.7(2.5) 

52.1(2.0) 

60.8(2.4) 

E5 

0.0926(15) 

51.9(2.1) 

59.6(2.4) 

53.6(2.4) 

62.4(2.7) 

F6 

0.0562(9) 

50.0(2.5) 

58.1(2.7) 

49.6(2.0) 

57.4(2.3) 

F7 

0.0449(8) 

46.1(1.9) 

54.1(2.2) 

47.9(2.0) 

56.1(2.2) 

G8 

0.0260(5) 

48.1(2.9) 

53.7(2.6) 

- 

- 

N5 

0.0940(19) 

52.6(2.8) 

57.7(3.6) 

- 

- 

N6 

0.0662(10) 

51.1(3.0) 

53.3(2.9) 

51.2(2.4) 

55.8(2.4) 

07 

0.0447(7) 

51.7(3.7) 

54.4(3.1) 

- 

53.7(2.4) 

LO-a2 

y®’'P,a = 0 

45.2(2.9) 

43.6(2.4) 

LO-ai 

= 0 

42.0(4.2) 

37.8(3.3) 

NLO-a^ 

y®’^P,a = 0 

44.9(2.9) 



NLO-a^ 

2 /®’=P,a = 0 

41.7(4.2) 




Table 4: Hyperfine splittings in the light- and strange-quark sector. 


For the hyperfine spin splittings, we obtain 

me* - me = 41.7(4.2)(3.2)a(0.3)x MeV , (3.3a) 

m-B* - me, = 37.8(3.3)(5.8)^ MeV , (3.3b) 

where we quote the mean value obtained with Cd = 1 for the B-meson case and add 
the difference w.r.t. the result obtained from the ansatz with Cd = 0 as a systematic 
error estimate for the chiral extrapolation. As mentioned earlier, we also account here 
for a systematic error between linear and quadratic continuum extrapolations. The 0(a) 
extrapolations are shown for both the B and the Bg system in Fig. 3. The filled/empty 
symbols at the physical point are the results using either an 0(a), or an O(a^) term in the 
continuum extrapolation. 

Our result for the hyperfine splitting for the B-system is in good agreement with the 
experimental value mp* — rnp = 45.78(35) MeV [1], whereas the Bg hyperfine splitting 
differs noticeably from the experimental value mp* — m-Ba = 48.7_2\ MeV. Our results 
are smaller than the experimental value and than our result for the hyperfine splitting 
of the B (the opposite of the situation for the experimental values). Since the hyperfine 
splitting came out far too small in the quenched approximation [39], this is suggestive 
of a residual quenching effect from the quenching of the strange quark in our Af = 2 
simulations. Moreover, Fig. 3 indicates larger cutoff effects in the case of the hyperfine 
splitting for the Bg. 

3.2 Excited states 

For the mass splittings between the ground state B,j and the first excited state, denoted 
here by B'’ and Bg, we obtain 

m-Qr — m-Q = 791(73) MeV , [m-Qr — 7713 ]®*’^* = 701(65) MeV , (3.4a) 

rriBi — ruBs = 566(57) MeV , [mB| — = 547(34) MeV . (3.4b) 


11 










Ar.tjm 

[GeV] 


[GeV] 

e-id 

y 

HYPl 

HYP2 

HYPl 

HYP2 

A4 

0.0771(14) 

0.692(37) 

0.703(33) 

0.619(27) 

0.625(25) 

A5 

0.0624(13) 

0.736(18) 

0.738(17) 

0.655(16) 

0.664(15) 

B6 

0.0484(9) 

0.622(66) 

0.665(57) 

0.552(56) 

0.592(52) 

E5 

0.0926(15) 

0.753(27) 

0.758(25) 

0.676(24) 

0.682(23) 

F6 

0.0562(9) 

0.719(57) 

0.721(54) 

0.658(52) 

0.654(47) 

F7 

0.0449(8) 

0.760(32) 

0.756(26) 

0.678(32) 

0.676(30) 

G8 

0.0260(5) 

0.67(11) 

0.70(9) 

0.63(11) 

0.65(9) 

N5 

0.0940(19) 

0.83(11) 

0.71(14) 

0.697(35) 

0.709(31) 

N6 

0.0662(10) 

0.78(9) 

0.76(8) 

0.634(72) 

0.668(58) 

07 

0.0447(7) 

0.71(12) 

0.71(8) 

0.64(10) 

0.63(8) 

L0-a2 

y’^^P,a = 0 

0.787(71) 

0.701(65) 

LO-ai 

y^^P,a = 0 

0.84(12) 




Table 5: Mass gaps between the excited state BJj- and the Bj-meson at static and next-to-leading 
order HQET. 




Ar,sW 

[GeV] 


[GeV] 

e-id 

y 

HYPl 

HYP2 

HYPl 

HYP2 

A4 

0.0771(14) 

0.774(45) 

0.748(34) 

0.640(23) 

0.650(22) 

A5 

0.0624(13) 

0.760(15) 

0.753(15) 

0.670(11) 

0.677(11) 

B6 

0.0484(9) 

0.719(49) 

0.704(35) 

0.641(22) 

0.641(21) 

E5 

0.0926(15) 

0.791(66) 

0.676(44) 

0.609(29) 

0.624(25) 

F6 

0.0562(9) 

0.689(30) 

0.722(24) 

0.643(18) 

0.642(17) 

F7 

0.0449(8) 

0.665(24) 

0.660(19) 

0.593(15) 

0.601(14) 

N6 

0.0662(10) 

0.698(49) 

0.718(36) 

0.621(38) 

0.643(32) 

07 

0.0447(7) 

- 

0.736(38) 

- 

0.669(29) 

LO-a^ 

y^^P, 0 = 0 

0.570(47) 

0.547(34) 

LO-ai 

y^^P, 0 = 0 

0.519(74) 




Table 6: Mass gaps between the excited state Bg- and the Bg-meson at static and next-to-leading 
order HQET. 
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e-id 

y 

A l/m 
^r.d 

[MeV] 


[MeV] 

HYPl 

HYP2 

HYPl 

HYP2 

A4 

0.0771(14) 

73(27) 

78(21) 

133(38) 

98(26) 

A5 

0.0624(13) 

71(11) 

75( 8) 

90(11) 

76( 9) 

B6 

0.0484(9) 

70(31) 

73(20) 

78(45) 

64(28) 

E5 

0.0926(15) 

77(17) 

75( 9) 

182(71) 

52(40) 

F6 

0.0562(9) 

61(31) 

68 (21) 

47(24) 

80(17) 

F7 

0.0449(8) 

82(18) 

80(11) 

71(19) 

59(15) 

G8 

0.0260(5) 

36(26) 

54(15) 

- 

- 

N5 

0.0940(19) 

132(110) 

6(95) 

- 

- 

N6 

0.0662(10) 

149(59) 

94(31) 

77(31) 

76(18) 

07 

0.0447(7) 

75(25) 

84(17) 

- 

67(26) 

LO-ai 

y^^P,a = 0 

0.090(40) 

0.019(46) 


Table 7: Subleading HQET contributions to the mass gaps between excited and ground states. 




Figure 4: Chiral and continuum extrapolation of — niB (left) and mer — (right). Both 
represent an extrapolation with ansatz (2.21) and E”g = 0 for HQET to order 1/mb. The open 
triangle represents the corresponding result for extrapolating the static order data. For details see 
Tables 5 and 6. We also add the continuum point (filled black pentagon) corresponding to the 
combination given in eq. (3.4). 
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Figure 5: The measured gap mer — me, vs. 2mT^ for the excited states with q = d (left) 
and q = s (right) on each individual CLS ensemble and with HYP2 discretization of the static 
quark. The boxes are an estimate of the expected gap for a two-hadron state B*(—p) -|- 7r(p) 
with one unit of lattice momentum p = |p| = 271 jL. The lower edge of the boxes corresponds to 
-I- -I- (277 jVf- and the upper edge includes an additional contribution p^/(2tob*) as a 

naive estimate for the kinetic energy of the B*. The dotted line is A = 27X1^^ corresponding to a 
three-hadron state B^ -|- tt -|- tt. 


after adding the individually extrapolated results for [m,Br — and [mB^ — 

In Fig. 4 these values are shown (as pentagons, slightly shifted at the physical point) for 
comparison with an extrapolation according to eq. (2.21). The raw data is collected in 
Tables 5, 6, and 7. 

We conclude this Section by remarking that for the excited states the interpretation 
of our results in terms of mass differences of physical one-meson states, e.g. “radial exci¬ 
tations”, is not straight forward. Although our values for these mass gaps are larger than 
what a multi-hadron state made of, e.g., a B|*j-meson and a small number (< 2) of phys¬ 
ical pions would produce, we cannot unambiguously conclude that our states actually 
correspond to radial excitations of the ground-state Bg-mesons. In a rigorous approach, 
states above multi-hadron thresholds need to be treated as resonances in order to obtain 
precise values for their masses and widths, and to associate them with the states observable 
in experiments. 

Since our lattice study is unquenched, the excited states can also be multi-particle 
states involving additional pions, beside the desired one-particle state. While it has been 
argued that the overlap of single-hadron interpolating operators to multi-hadron states is 
small [46,47], the two-hadron states may have a weaker volume suppression [48]. Moreover, 
a chiral extrapolation linear in might not be adequate for a multi-hadron state, and a 
different extrapolation ansatz for these states, e.g. linear in might also yield somewhat 
smaller mass gaps. 

In continuum (and infinite-volume) physics, any excited Bg-meson state, B^, which 
strongly decays into an Awave Bg-bvr state, i.e. a resonance in the Awave Bg-bvr scattering 
channel, implies that the corresponding two-hadron state B^ -|- tt with relative angular 
momentum i has the correct quantum numbers to couple to our interpolating operators 
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Figure 6: The same data points as in the right panel of Fig. 5 for the measured mass gaps, 
together with estimates of the energy gap for a two-hadron state B*(—p) -|- K(p) (boxes) and a 
three-hadron state B^ -I- K -|- tt (dashed lines). 


used for the Bg. For the q = d sector, the set of possible two-hadron excited states includes 

• B* -|- vr in P-wave which would have a non-interacting two-hadron energy gap of 
A > 181 MeV, 

• B'* -|- vr in P-wave, where B'* is a radial excitation (= 1“), 

• Bg -|- TT in S-wave, where Bg is an orbital excitation (J^ = 0^), 

• B 2 -|- vr in D-wave, where B 2 is an orbital excitation (J^ = 2'*'), like the observed 
B2(5747)‘^ state which would lead to a non-interacting two-hadron energy gap of 
A > 598 MeV. 

On a finite lattice, the energies of states with non-zero relative angular momentum are 
lifted due to the minimal momentum of |p| = 27r/L. Since the s-quark is quenched in our 
study, we cannot have excited states in the B sector with an s-quark from the sea, like 
B* -|- K. On the other hand, in the Bg sector there are two flavour combinations for each 
of the corresponding excited state of the B sector, e.g. B* -|- vr and B* -|- K (corresponding 
to B* -|- vr). 

In Figs. 5 and 6 we show our data points on the various ensembles together with an 
estimate for the gaps of some two-hadron states. The comparison indicates that some of 
the two-hadron states can be close in energy to the excited states we measured. Depending 
on the pion mass on a given ensemble, the energy gaps determined according to eqs. (2.16) 
may in fact be the energy splitting to the lowest lying B^ -|- vr state, if this is lighter than 
the radial excitations B^. 

Below three- or more-hadron thresholds, the infinite volume scattering matrix (up 
to corrections which vanish exponentially in the spatial lattice extent) may be inferred 
from the finite volume energy spectrum [49-53]. In practice this requires the construction 
of correlation matrices which contain two-hadron as well as single hadron operators. To 
date, this procedure has been carried out successfully in some simple systems such as vr—vr 
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Source 

me. 

me- - tob 

ms* - 

wb, - mB 

niBr - tub 

rUBl ~ 

i 

5383(63) 

41.7(4.2) 

37.8(3.3) 

88.9(5.7) 

791(73) 

566(57) 

A4 

0.46 % 

2.52 % 

3.51 % 

12.48 % 

9.79 % 

9.93 % 

A5 

0.30 % 

1.32 % 

1.81 % 

3.14 % 

6.67 % 

5.78 % 

B6 

0.03 % 

1.50 % 

0.39 % 

0.23 % 

1.04 % 

0.06 % 

E5 

0.28 % 

0.40 % 

1.80 % 

0.46 % 

4.72 % 

1.51 % 

F6 

0.10 % 

0.34 % 

0.45 % 

8.62 % 

5.07 % 

6.72 % 

F7 

0.21 % 

0.82 % 

1.50 % 

34.56 % 

17.87 % 

23.42 % 

G8 

0.53 % 

5.25 % 

0.00 % 

0.05 % 

6.46 % 

0.03 % 

N5 

1.90 % 

8.04 % 

0.45 % 

0.01 % 

7.94 % 

0.05 % 

N6 

5.97 % 

1.33 % 

29.18 % 

31.56 % 

14.60 % 

27.70 % 

07 

4.50 % 

6.34 % 

13.49 % 

8.05 % 

25.32 % 

14.42 % 

UJ 

62.84 % 

31.75 % 

46.61 % 

0.04 % 

0.09 % 

0.05 % 

Za 

21.13 % 

0.27 % 

0.63 % 

0.75 % 

0.40 % 

0.19 % 


Table 8: Distributions of relative squared errors, af / among different sources i for the mean 

values given in (3.1)-(3.4). In the second row we repeat them in MeV for the reader’s convenience, 
including only the statistical error of the mean. Remaining relative errors enter through the scale 
setting procedure only. 


and vr—K scattering, and the energy dependence of the scattering phase shift has been 
determined with a snfhcient resolntion to clearly discern resonant behavionr (see e.g. [54] 
for a recent review). Withont performing snch a dedicated stndy, we are not able to 
conclnsively determine the natnre of the excited state we observe. 

4 Conclusions 

In this paper, we have presented resnlts for the B-meson spectrnm obtained in the frame¬ 
work of lattice HQET expanded to 0(l/mb). Within this approach the existence of a 
continnnm limit is gnaranteed, as nnmerically tested with high accnracy in previons stnd- 
ies [55,56]. In contrast to the HQET expansion in continnnm pertnrbation theory, onr 
approach is manifestly non-pertnrbative in the strong conpling. We perform the contin¬ 
nnm extrapolation from lattice resolntions in the range 0.08-0.05 fm. Pion masses, in onr 
setnp with Nf = 2 degenerate flavonrs, reach down to valnes of abont 190 MeV. The accn¬ 
racy is at the 10% level, for the different splittings presented (e.g., hyperhne and Bj-Bs 
splittings), and we always hnd consistency within two standard deviations with valnes from 
the PDG, whenever a comparison is possible. 

Hyperhne splittings probe higher-order terms in HQET and the reported resnlts rep¬ 
resent an important check on the validity and the reliability of the asymptotic HQET 
expansion, trnncated at NLO, at the b-qnark mass scale. Compared to previons qnenched 
resnlts, we observe a signihcant shift for the hyperhne splitting in the B-meson sector, 
which now agrees with the experimental determination. In onr Nf = 2 simnlations, the 
hyperhne splitting in the Bg sector appears to snffer from a residnal qnenching effect, which 
is in line with what was seen for N{ = 0. In order to ascertain that the qnenching of the 
strange qnark is indeed the root canse of the rednced Bg hyperhne splitting seen here, we 
plan to extend onr compntations to simnlations of the N{ = 2 + 1 theory [57]. 

A dominant sonrce of nncertainty in onr resnlts is represented by cntoff effects (see 
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Tab. 8). This does not come unexpected since we have not implemented 0(a) improvement 
at 0(l/mb). While implementing a fully non-perturbative improvement programme at 
this order is probably too difficult, one may consider perturbative (tree level or one-loop) 
improvement for future applications. 

We also determined the mass gaps for excited states in both the B and the Bg sectors. 
The results are consistent with a radial splitting, e.g. as computed for the Be system in [9], 
but these excited states might also be two- or multi-hadron states. 

Knowledge of the mass splittings is relevant for the computation of hadronic pa¬ 
rameters within the sum-rules approach and when comparing results from the lattice to 
sum-rules estimates. The mass gaps of excited states are also an important information 
for the computation of form-factors on the lattice, for example for the B —)• 'K^v and the 
Bg —)• K£i/ decays, as currently endeavored by the ALPHA Collaboration [12,14], and in 
general in the spectral analysis of two- and three-point functions. 
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A Determination of the statistical error in the presence of autocorrela¬ 
tions 

In order to compute the statistical error in the presence of both correlations between the 
different observables and of autocorrelations along the HMC trajectories by which our 
ensembles were generated, we employ the methods of [58-61], which we briefly outline 
below. 

The starting point is the computation of the “primary” observables C/,, where i labels 
the Ameas gauge conhgurations, and a is an aggregate label for the different correlators 
measured (stat, spin, kin), the Euclidean time separation t between the source and sink, 
and the different smearing levels employed at source and sink. The gauge average Ca and 
the variance are computed as usual. To estimate the true statistical error of the 
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gauge average, we also require the integrated autocorrelation time which is computed 
from the autocorrelation function 

rS ^ rc.c,(T) = if; (cj- - c„) (c‘ - Cj), (A,i) 

i=l 

where r is the separation in simulation time along the Markov chain. In order to take into 
account the long-time tail of F, the conservative estimate [59] 


..int _ 1 1 

- 2 + 


ri^i(o) 



(A.2) 


is used, where is an estimate of the exponential autocorrelation time of the Markov 
chain. The values used in our analysis are listed in Table 2. The window size W is 
automatically chosen as the point t = W where rQo(r) comes close to zero within about 
1.5 of its estimated error. The true statistical error is then given by 



2r, 


int Oa 
Co 




(A.3) 


For derived observables which are functions of gauge averages of the primary 
observables Cq, and in our case include the generalized eigenvalues and eigenvectors as 
well as the energies derived from them, we compute the derivatives 


j = ^ 

dC^ 


and the autocorrelation function [59] 


o,/3 


(A.4) 


(A.5) 


The variance of the derived observable Do,/ is then given by 



_ p{2) 

Ol'OL 


.(0) 


(A.6) 


and its statistical error by 


2 o.^int ^ 


D' 


(A.7) 


where the integrated autocorrelation time tV}} is again estimated using (A.2), with T^^^ 
substituted for T^^^. 

Since the extraction of plateau values from a weighted fit requires knowledge of the 
errors of the individual points, in principle this procedure should be iterated, with the 
plateau averages as secondary derived observables of the derived observables. However, 
the integrated autocorrelation times for the effective energies in the plateau region do not 
markedly differ. Therefor, it is sufficient to employ their variances (as estimated using 
e.g. a Jackknife procedure for error propagation) to weight the fit, treating only the fitted 
values as derived observables. This simplified procedure has been adopted here. 

In order to extract the final answer in physical units, the plateau values must be com¬ 
bined with each other and with the HQET parameters and lattice spacing, propagating 
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the errors on each of those to the final resnlt, where the HQET parameters are statisti¬ 
cally independent of the large-volnme observables; similarly, in extrapolating to the chiral 
and continnnm limits, the resnlts obtained from different ensembles are statistically in¬ 
dependent, and their contribntion to the error of the final resnlt / can thns be added in 
qnadratnre: 




E 




^ df ^ df 

dYj ’ 


(A.8) 


where the Yi are the additional parameters, Cy.y^. is their (known) covariance matrix, and 
(Tj(e) is the error compnted according to eqn. (A. 7) when taking into acconnt only the 
flnctnations of Ca on ensemble e [59-61]. 


B Effective energies and matrix elements 


In this Appendix we provide the nnmerical resnlts of onr GEVP analysis after performing 
a weighted plateanx average 




w{t) = (B.l) 


for t G [t min ■ t mav ] and to > tmin/2, see Section 2.2. The specific valne of tmax is irrelevant 
dne to noise dominated data at large time separations.^ As a conseqnence of the expo¬ 
nential growth of the noise-to-signal ratio, the qnoted errors are dominated by the error 
at tmin and we decided to qnote tmin as snbscript to the statistical error in the following 
Tables. For determining tmiiu we nse to = ^ ~ 1 ia the GEVP, and then tg = tmin ~ 1 once 
tmin is fixed. The errors qnoted are those entering the corresponding effective energy plots 
presented below and do not contain the tail contribntion of the error. Since the antocorre- 
lation is expected to be the same among different time slices in the platean region, it does 
not affect the estimate (B.l) and can be added, inclnding the exponential tail, whenever 
needed explicitly. In fact, all valnes qnoted for derived observables as presented in the 
main text have the exponential tail inclnded as discnssed in Appendix A. 

The same procednre has been nsed to obtain matrix elements p^. We inclnde them in 
the present paper (Table 13-14) in order to complete the data set nsed in [15,41]. 

A special case is the ensemble labeled A5 for which we have two independent Monte 
Garlo histories/replicas (A5c and A5d). Here, measnrements and a snbseqnent GEVP 
analysis have been independently performed for both replicas in the heavy-light sector 
bnt only on A5d for heavy-strange. The resnlts from different replicas are then combined 
before derived observables as presented in the main text are compnted. 
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Figure 7: Effective energies following our GEVP analysis in the heavy-light (top) and heavy- 

strange (bottom) sector on ensemble A4. 
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Figure 8: Effective energies following our GEVP analysis in the heavy-light (top) and heavy- 

strange (bottom) sector on ensemble A5c. 
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Figure 9: Effective energies following our GEVP 

ensemble A5d. 
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Figure 10: Effective energies following our GEVP analysis in the heavy-light (top) and 

heavy-strange (bottom) sector on ensemble B6. 
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Figure 11: Effective energies following our GEVP analysis in the heavy-light (top) and 

heavy-strange (bottom) sector on ensemble E5g. 
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Figure 12: Effective energies following our GEVP analysis in the heavy-light (top) and 

heavy-strange (bottom) sector on ensemble F6. 
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Figure 13: Effective energies following GEVP analysis in the heavy-light (top) and heavy- 

strange (bottom) sector on ensemble F7. 
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Figure 14: 

ensemble G8. 



Figure 15: Effective energies following onr GEVP analysis in the heavy-light sector on 

ensemble N5. 
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Figure 16: Effective energies following our GEVP analysis in the heavy-light (top) and 

heavy-strange (bottom) sector on ensemble N6. 
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Figure 17: Effective energies following our GEVP analysis for heavy-light (top) and heavy- 

strange (bottom) sector on ensemble 07c. 
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